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Abstract. The TeV emission of low power BL Lac objects has been established by the detection of an handful of them. 
The knowledge of the level of the TeV emission and its spectrum can shed light on the particle acceleration mechanisms, 
and it is especially important to assess the still uncertain level of the far infrared background radiation, which can absorb the 
TeV photons through photon-photon interactions. In view of these implications, it is necessary to enlarge the number of TeV 
detected sources, and to find them at different redshifts. To this aim, we propose a general and simple criterium to select the 
best TeV candidates, and produce a list of them with flux estimates above 40 GeV, 300 GeV and 1 TeV. 
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1. Introduction 

The discovery that blazars (Flat Spectrum Radio Quasars and 
BL Lac objects) are very strong 7 -ray emitters have renewed 
the interest about them, and opened new perspectives in the 
comprehension of the physics of these objects. 

The observations by EGRET (Hartmann et al. 1999), on¬ 
board the Compton Gamma Ray Observatory, led to the discov¬ 
ery that blazars emit most of their power in the 7 -ray band, and 
that their Spectral Energy Distribution (SED) is characterized 
by two broad peaks, now commonly (but not unanimously), 
interpreted as due to synchrotron and inverse Compton (IC) ra¬ 
diation, respectively (e.g. Maraschi et al. 1992; Dermer et al. 
1992; Sikora et al. 1994; Ghisellini & Madau 1996; but see 
Mannheim 1993; Rachen 1999; Muecke & Protheroe 2000; 
Aharonian 2000 for a different interpretation). Eor the first time 
we were able to study their entire SED in a comprehensive 
way, finding differences among subclasses of blazars about the 
frequency location of the two broad peaks, their relative lumi¬ 
nosity and variability behaviors in different bands. Considering 
the EGRET sources and three complete blazar samples, it was 
found a correlation between the location of the two broad peaks 
and the observed bolometric luminosity (Eossati et al. 1998, 
hereafter F98; Ghisellini et al. 1998). Blazars seem to form a 
well defined sequence, with low powerful objects having both 
peaks at a similar level of luminosity, and located at higher fre¬ 
quencies than in more powerful objects, in which the IC peak 
dominates the emission. In the BL Lac class, the hrst kind of 
sources were named High frequency Peak BL Lacs (or HBL 
for short) by Padovani & Giommi (1995), while the latter sub- 
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class was called Low frequency Peaked BL Lacs (LBL). Recent 
observations of high redshift {z > 4) blazars (Eabian et al. 
2001a, 2001b; Celotti 2001), and of low power BL Lac objects 
(Costamante et al. 2001a, 2001b) have extended the blazar se¬ 
quence at both ends, resulting in agreement with the original 
trend. 

At TeV energies, the detection and study of blazar objects 
by ground based Cherenkov telescopes have been limited up 
to now to few HBL sources (Mkn 501, Mkn 421, PKS 2155- 
304, lES 2344H-514, see Catanese & Weekes 1999), though 
disclosing new and fundamental aspects of the blazar behav¬ 
ior. These observations monitor the behavior of the most en¬ 
ergetic electrons of the source, thus shedding light on the ac¬ 
celeration mechanism working at the most extreme conditions. 
The very rapid variability observed at these energies (Mkn 421 
doubled its TeV flux in less than 20 minutes, see Gaidos et 
al. 1996, Catanese & Weekes 1999), coupled with the require¬ 
ment that the source must be transparent with respect to the 
photon-photon process, tightly constrains the physical parame¬ 
ters, such as the source size and its beaming factor. In addition, 
Mkn 421 showed a tight correlation between the emission in 
the X-ray and TeV bands (Maraschi et al. 1999, Takahashi et 
al. 2000, Krawczynski et al. 2001), implying that the radiation 
produced in the two bands is co-spatial and produced by the 
same electrons; this is of crucial importance to constrain any 
emission model. 

The strong connection between the TeV and X-ray emis¬ 
sion was also clearly evident during the 1997 flare of Mkn 
501, when this source was observed by the X-ray satellite 
BeppoSAX in an extreme spectral state, with a synchrotron 
peak frequency close to 100 keV or even more (Plan et al. 
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1998). Mkn 501 was found to have increased at least tenfold 
its luminosity, with most of it radiated at high X-ray energies. 
At the same time, the source underwent a major flare in the TeV 
band (Catanese et al. 1997b; Aharonian et al. 1997; Protheroe 
et al. 1997; Djannati-Atai et al. 1999), and continued to be ac¬ 
tive (and well visible) in the TeV band for several months. This 
dramatic behavior can be explained by a synchrotron inverse 
Compton model, taking into account the effects introduced by 
the Klein Nishina scattering cross section and the constraints 
posed by the transparency of the source with respect to photon- 
photon collisions producing electron-positron pairs. It is not 
clear if the simultaneous variations in the X-ray and TeV bands 
can be completely accounted for by a simple one-zone homo¬ 
geneous synchrotron self-Compton model (see Tavecchio et al. 
2001 ), or if we need some extra and more quiescent source of 
IR-optical photons (i.e. Ghisellini 1999). 

TeV observations of BL Lacs are also particularly inter¬ 
esting because, being the only extragalactic sources known to 
emit at these energies, allow an independent estimate of the ex¬ 
tragalactic IR background (IRB). Direct measurements of the 
IRB are affected (up to now) by relatively large uncertainties, 
due to the heavy contamination of foreground objects (for a 
review see Hauser & Dwek 2001). Because TeV photons can 
be absorbed by IR photons for the pair production mechanism, 
the analysis/study of high energy spectra from sources at dif¬ 
ferent redshifts gives an independent measure of the IRB level 
(Stecker et al. 1992). 

These kind of studies are just started, and the first conclu¬ 
sions are based on the two most observed TeV BL Lacs, namely 
Mkn 421 and Mkn 501. The main uncertainty here is the knowl¬ 
edge of the primary blazar spectrum, which could have an in¬ 
trinsic cut-off at high energies, or could be affected by absorp¬ 
tion due to IR photons produced locally. A first and preliminary 
confirmation of the IRB absorption has been obtained compar¬ 
ing the spectra of Mkn 421 and Mkn 501, showing a cut-off 
in the power-law spectra at approximately the same energies 
(Krennrich et al. 2001), as expected since they have similar red- 
shifts. 

The direct measurements of the IRB flux lead to pre¬ 
dict quite a strong absorption at TeV energies. If true, this 
in turn would imply an unusual primary spectrum above ~ 
1 TeV, which must have an excess above the extrapolation 
from lower energies, leading to the so-called IRB-TeV puz¬ 
zle (see e.g. Protheroe & Meyer 2000; Aharonian, Timokhin & 
Plyasheshnikov, 2001; Berezinsky 2001) []. 

However, to draw unambiguous conclusions, we need an 
ensemble of sources located at different redshifts and a detailed 
knowledge of the X-ray flux and spectrum; being produced by 
the same electrons (in synchrotron inverse Compton models), 
this would help in predicting the shape of the TeV emission. 

To this aim we consider in this paper several published 
samples of bright BL Lac objects, for a total of 246 different 

* It has been proposed that this “puzzle” could be solved by 
quantum-gravity theories predicting the breaking of Lorentz invari¬ 
ance: as a consequence there should be a modification in the energy 
threshold for the 7-7 —> process, explaining why TeV photons are 
not heavily absorbed by the IRB (Amelino-Camelia & Piran 2001). 


objects, and propose a simple and handy tool to identify and 
select the most promising candidates. The main point we em¬ 
phasize concerns the requirement of both high energy electrons 
and sufficient seed photons to originate the TeV emission. We 
therefore consider as best candidates those BL Lac objects hav¬ 
ing not only their synchrotron peak located at high energies, but 
also having sufficient radio-through-optical flux. 

We therefore expand the work first done by Stecker et al. 
(1996), concerning only the Einstein Slew survey sample of 
BL Lac objects, both by considering other samples and also by 
introducing a different, albeit still simple, criterium to identify 
the best candidates. 

2. BL Lac objects detected at TeV energies 

At present, only few blazars have been detected by Cherenkov 
telescopes. In Table we report the basic properties, fluxes 
and detection thresholds for those objects detected at least 
once above 5 (t (see Catanese & Weekes 1999 and references 
therein). We have also added lES 1426H-428, whose detection 
was reported at about the 4.5a level by WHIPPLE (Horan D. 
et al., 2000). This was confirmed by CAT (Djannati-Atai, priv. 
comm.), by WHIPPLE again in 2001 (Horan D. et al., 2001), 
and by HEGRA (Aharonian 2001b). The latter two detections 
were at the 5cr level. 

Of these objects however, only Mkn 421 and Mkn 501 
have been firmly confirmed by repeated detections by many 
Cherenkov telescopes, and since 1995 they have been exten¬ 
sively monitored. Thanks to their relatively high fluxes and re¬ 
peated activity states, these are also the only two sources with 
spectral informations up to ~ 17 and 20 TeV (Krennrich et al. 
2001 ). 

Two other objects have been reported to be TeV emitters, 
but with a somewhat lower significance. lES 1959H-650 was 
observed in 1998 by the UTAH Seven Telescope Array, and 
found with an excess of 3.9 (t, reaching in two periods ~ 5cr, but 
no fluxes were reported (Nishiyama et al. 1998). 3C 66A was 
instead detected once at ^ 4.2cr by the Crimean Astrophysical 
Observatory in 1996 (Neshpor et al. 1998), with a detection 
threshold at 0.9 TeV. If confirmed, this source could be ex¬ 
tremely important, since due to the high redshift (z=0.444), its 
TeV emission should be suppressed according to the actual es¬ 
timates on the 7-7 ^ optical depth due to the IRB (Stecker 
et al. 1999). 

3. The Synchrotron self Compton process and 
TeV emission 

The BL Lac objects so far detected at TeV energies are rela¬ 
tively low powerful blazars, with no broad emission lines and 
no signs of thermal emission (i.e. the blue bump) produced by 
the accretion disk. This suggests and favors a synchrotron self- 
Compton (SSC) origin of the TeV flux. However this may be an 
over-simplification, since, besides the jet region containing the 
energetic electrons responsible for the high energy emission, 
other sites, both in the jet and externally to it (e.g. a molecu¬ 
lar torus, a thin scattering plasma surrounding the jet, or the 
walls of the jet itself), may be important in producing the soft 
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Fig. 1. Spectral energy distributions of the 5 BL Lac objects already detected at TeV energies. The solid and dashed lines refer to 
the SSC model described in Section 3.2 and to the SED constructed using the parameterization described in Fossati et al. 1998, 
with the modification described in this paper. The input parameters used for the SSC models are given in Ghisellini, Celotti & 
Costamante (2001). 


Source 


Tsghz 

Jy 

T 5500 A 

mjy 

ClkeV 

+Jy 

-^TeV 

Eth 

TeV 

Signif. 

Tel. 

Mkn421 

0.031 

0.722 

9.75 

9.98 

1-80 

0.3 

conf. 

WHIPPLE, CAT, HEGRA 

lES 1426+428 

0.129 

0.038 

1.56 

6.83 

0.9 

0.4 

4.5 — 5(7 

WHIPPLE, (CAT), HEGRA 

Mkn 501 

0.034 

1.371 

10.31 

9.44 

0.8-49 

0.3 

conf. 

WHIPPLE, CAT, HEGRA 

PKS 2155-304 

0.116 

0.310 

22.34 

14.6 

4.2 

0.3 

6.8(t 

Durham Mark6 

lES 2344+514 

0.044 

0.215 

3.54 

2.91 

1.7 

0.35 

5.2(7 

WHIPPLE 


Table 1. Data for the BL Lac objects detected in the TeV band. Significances are listed only for unconfirmed sources. Data from 
Catanese & Weekes (1999) and Krennrich et al. (2001), and references therein. 1) integrated flux in units of 10“^^ photons cm“^ 
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seed photons to be scattered at high energies. However, a sim¬ 
ple one-zone and homogeneous SSC model can account for 
the observed SED of all the TeV BL Lacs, and due to its sim¬ 
plicity it is the basic framework we will use in the following. It 
is also the model assumed by the previous work on this issue, 
namely the paper by Stecker et al. (1996), which we discuss 
below, before introducing our views. 

3.1. The Stecker et al. (1996) scenario 

Based on the SSC framework, Stecker, De Jager & Salamon 
(1996) have used simple scaling arguments to predict the 7 - 
ray fluxes for HBL objects, and so to select good candidates 
for TeV emission. According to the SSC model, the inverse 
Compton component has a spectrum which is similar to the 
synchrotron one (both roughly parabolic on a logarithmic i/F^, 
plot, Macomb et al. 1995, L98), but upshifted by ~ 7peak 
the Thomson regime), where 7 peak is the Lorentz factor of the 
electrons emitting at the peak. Using as “template” the SED 
of Mkn 421, at the time the only HBL source detected both 
at GeV (EGRET) and TeV (WHIPPLE) energies (i.e. with in¬ 
formations on both sides of the Compton peak), they found an 
upshifting factor of ^ 10®, and a Lq /Esyn ~ 1- Assuming then 
for simplicity that all HBL objects have the same properties as 
those found for Mkn 421, the found upshifting factor allowed 
to derive the following scaling law: 

t'o^o t^GevPCeV , J^xFx I^TevFTeV ,,, 

7 - -- 7 - and ---- - -, ( 1 ) 

r^syn xjC J^syn J^C 

where Fo is the monochromatic flux at the optical frequency 
i/o- Given Lc/^syn ~ 1, a direct relation for the energy fluxes 
is obtained, in the GeV and TeV ranges: 

’^GevFcgy ^ VoFo and J^TeV.PreV ~ '^xFx ( 2 ) 

Quantitative estimates on the integral fluxes were then made 
using the average spectral index for BL Lacs between 0.1 and 
10 GeV {a = 0.8), and the Mkn 421 slope above 0.3 TeV 
(a = 1 . 2 ). 

In this scenario, therefore, the main selection criterium is 
the X-ray flux level: the energy flux around 1-2 keV gives 
directly the expected level of the emission roughly around 0 . 2 - 
0.3 TeV. In this context, sources with the highest synchrotron 
t'peak are expected to emit more in the TeV band, since the 
Compton peak progressively moves into the TeV window. 

3.2. Energetic electrons and seed photons 

The two ingredients for the formation of a strong TeV emis¬ 
sion by the synchrotron self Compton process are the density 
of electrons energetic enough to produce TeV photons and the 
density of seed photons to be scattered. Lor not too large syn¬ 
chrotron peak frequencies, the scattering process between pho¬ 
tons at the synchrotron peak and the electron producing the 
peak itself is in the Thomson regime. In this case the relevant 
electrons can scatter the bulk of the synchrotron photons to 
high energies. 



Fig. 2. SEDs calculated with a simple one-zone SSC model as 
described in Section 3.2, in which 7 peak increases (as labeled). 
All other parameters are kept constant. 

For larger synchrotron peak frequencies (and assuming the 
same value of the magnetic held) we have more energetics elec¬ 
trons in principle capable to produce photons of larger energies. 
However in this case the Klein-Nishina decline of the scatter¬ 
ing cross section is important, disfavoring those scatterings be¬ 
tween the relevant electrons and the synchrotron photons ly¬ 
ing at the peak. In this case the scattering process producing 
TeV photons is produced by seed photons of lower energies, 
with a corresponding reduced synchrotron energy density. In 
other words, not all the synchrotron photons are used to form 
the Compton spectrum, but only the photons of energy (in the 
comoving, primed, frame) hF < rrieF can efficiently (i.e. 
they scatter in the Thomson regime) contribute to the Compton 
emission. There is then a trade-off: to produce large TeV fluxes 
we need very energetic electrons, corresponding to very large 
t'peak. but we also need enough seed photons, and this requires 
not extreme values of i^peak- 

To better illustrate this case we show in Fig. 2 a sequence of 
spectra derived by a standard SSC model in which all parame¬ 
ters but 7 peak are kept constant. In this model it is assumed that 
relativistic electrons between 71 and 72 are injected at a con¬ 
stant rate Qij) oc 7 “®-® [cm“® s“^], throughout a spherical 
source of radius R, magnetic held B and beaming factor 8. The 
steady state particle distribution is found through the continu¬ 
ity equation, accounting for synchrotron and inverse Compton 
radiative losses, electron positron pair production and Klein 
Nishina effects. Particles are assumed not to escape the source. 
In this case 7 peak is always coincident with 71 (see Ghisellini 
et al. 1998 for more details of the model). 

In Fig. 3 we report the flux integrated above a given fre¬ 
quency as a function of i^peak. corresponding to the spectra 
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Fig. 3. Top panel shows the flux integrated above the indicated 
energies corresponding to the spectra of Fig. 2. i^peak is the 
synchrotron peak frequency of the SEDs in Fig. 2. The bot¬ 
tom panel shows the radiation energy density available for 
scattering in the Thomson regime with electrons of energy 

TpeakTOeC^- 


shown in Fig. 2. In the bottom panel we also show the comov¬ 
ing radiation energy density which is available for scattering in 
the Thomson limit. From both Fig. 2 and Fig. 3 we can see that 
the level of the TeV emission initially increases for increasing 
7 peak, since in this case the scattering of the bulk of the syn¬ 
chrotron photons occurs in the Thomson regime, and the in¬ 
crease in 7peak results in more photons reaching TeV energies. 
This trend is reversed when > TOeC^/7peak, e-g- when 

r 

i^peak > 3.8 X Hz (3) 

1 - 1 - 2 : 

where i/peak = (4.3)j/L7peak'^/(l + '2^)’ I/L = eB/ (27rmec) 
is the Larmor frequency. The increase in the synchrotron lumi¬ 
nosity for increasing r'peak is due to the corresponding decrease 
in the self Compton component: since all the injected power is 
assumed to be radiated, the sum of the synchrotron and the self 
Compton luminosities must be constant. 

We conclude that very extreme BL Lacs, with very large 
values of t'peak. are not necessarily the best candidates to be 
strong TeV emitters since they can scatter relatively fewer pho¬ 
tons. The best TeV BL Lac candidates should be the one with 
both a large r'peak and a sufficiently strong soft seed IR photon 
emission. 


4. TeV candidate BL Lacs 


4.1. The Samples 

In order to select new candidates for TeV emission, we have as¬ 
sembled a catalogue of BL Lac objects using several published 
BL Lac samples, for which informations in all the three energy 
bands (radio, optical and X-ray) were available. We considered 
the Slew Survey Sample (Perlman et al. 1996), the Einstein 
Medium Sensitivity Survey (EMSS, Rector et al. 2000), the 
ROSAT All Sky Survey BL Lac sample (RASS, Bade et al. 
1998), the ROSAT All Sky Survey - Green Bank sample (RGB, 
Laurent-Muehleisen et al. 1999), the EXOSAT archive BL Lac 
catalogue (Giommi et al. 1990, Sambruna et al. 1994) and the 
1 Jy BL Lac sample (Shekel et al. 1993, Urry et al. 1996). We 
have also added all the HBL objects in Donato et al. (2001), 
who gives a list of all known blazars detected in X-rays for 
which a measure of the spectral index was available. 

The sources listed in more than one sample have been 
considered only once, and attributed with the following or¬ 
der: Slew, EMSS, RASS, RGB, EXOSAT, 1 Jy samples and 
those in the Donato et al. compilation (i.e. the EMSS sources 
in Fig. ^ and ^ are those not already in the Slew sample, the 
RASS sources are those not included in the EMSS and Slew 
samples, and so on). In this way we obtained an ensemble of 
246 different objects. The flux informations used are those re¬ 
ported in the respective catalogue papers (see each reference). 
Because the data were often provided in different formats, we 
have uniformed all the fluxes to monochromatic fluxes at 5 
GHz, 5500A and 1 keV, using the same spectral information 
adopted in the compilation of each sample, for consistency. In 
the radio band, the RASS sample reported the fluxes at 1.4 
GHz along with the radio spectral index, so we have calcu¬ 
lated the flux at 5 GHz assuming a power-law spectrum and 
the reported q;r. In the X-ray band we calculated the fluxes at 
1 keV from integrated fluxes (RASS and RGB) or fluxes at 2 
keV (Slew, EMSS) using the X-ray spectral indices there re¬ 
ported. The optical fluxes were calculated from the tabulated 
magnitudes, dereddened with the Ab values obtained from the 
NED database (Burstein & Heiles 1982). When the optical flux 
was not at 5500A (V filter), we obtained it extrapolating from 
the flux at the different effective wavelength, assuming a power 
law spectrum of slope aopt = 1- The NED database was also 
used to check the redshift for all sources, 65 of which do not 
have a reported value. 

For Fig. ^ and ^ the fluxes have been K-corrected with 
the respective catalogue spectral indices, assuming power-law 
spectra. When the slope was not reported, we used aR=o and 
ctopt = 1 for the K-correction of the radio and optical fluxes, 
respectively. The X-ray flux was K-corrected using the re¬ 
ported spectral indices. For the K-correction in case of un¬ 
known redshift we used z = 0.2. 

4.2. Selection of TeV candidate BL Lacs 

Fig. 4 shows the X-ray flux as a function of the radio flux for 
the BL Lacs objects in our ensemble. Note the locations of 
the already TeV-detected sources: they are among the brightest 
sources in both bands. This is not so obvious as it may seem 
at first sight, since a large X-ray to radio flux (hence, a lower 
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Fig. 4. BL Lac objects in the radio (5 GHz) and X-ray (1 keV) vF(y) plane. Sources belonging to different samples have different 
symbols, as labeled. The objects marked with hlled circles are those already detected at TeV energies (from left to right, lES 
1426+428, lES 2344+514, PKS 2155-304, Mkn 421, Mkn 501). Note that for these sources we have plotted two different X-ray 
states, connected by the vertical segment. The dotted lines delimiting the rectangle are at = 1.46/iJy and Tr = 31.6 mJy. 


radio emission for a given X-ray flux) indicates a large syn¬ 
chrotron peak frequency (see e.g. Fig. 8 in F98), which is the 
hrst requirement to emit in the TeV band (we must have many 
electrons energetic enough to emit copiously at TeV energies). 
Consider also that the radio emission (at 5 GHz) must be pro¬ 
duced in a large region of the jet (not to be self-absorbed), 
much larger than the part of the jet emitting high frequency 
radiation (as required by the very rapid variability). Therefore 
the link between the radio/X-ray emission and the TeV flux is 
more subtle. 

We have interpreted this property in the following way: 
to produce a large TeV flux by the IC process we need 


many electrons of random Lorentz factors 7 - 10^-10®. 
These electrons emit synchrotron photons of energies hv = 
1 .5i3(7/10®)^((5/10) keV, with B in Gauss. The seed photons 
most effective to interact with these electrons to produce TeV 
photons by the IC process are in the IR-optical band, since pho¬ 
tons of higher frequencies scatter in the Klein Nishina regime. 
We therefore propose that the radio flux measures the level of 
the relevant seed photons. If this is the case, then, for a given 
X-ray flux, sources that are brighter in the radio band are more 
likely to be TeV emitters. And indeed the objects already de¬ 
tected in the TeV band are bright both in the radio and in the 
X-ray bands. We think this is exactly the effect of the “trade- 
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Fig. 5. BL Lac objects in the optical (5500 A) and X-ray (1 keV) vFiv) plane. Sources belonging to different samples have 
different symbols, as labeled. The objects marked with hlled circles are those already detected at TeV energies (from left to 
right, lES 1426+428, lES 2344+514, Mkn 421, Mkn 501 and PKS 2155-304). Note that for these sources we have plotted two 
different X-ray states, connected by the vertical segment. The dotted lines delimiting the rectangle are at = 1.46/rJy and 
Fo = 0.49 mJy. The objects marked also with a black diamond are those inside the rectangle of the i'xFx-vb.Fr plot (Eig. 


off’ between high 7peak and the energy density of seed photons 
(see Sect. 3.2), since the objects with higher t'peak are, on av¬ 
erage, those with fainter radio emission for a given X-ray flux, 
i.e. those with lower (see the correlation between a^x and 
^'peak in F98, Wolter et al. 1998, Costamante et al. 2001). 

As shown in Fig. ^ in the region of high radio and X- 
ray flux around the detected TeV objects there are also other 
sources, that we therefore consider good candidates for TeV 
emission. The extension of such region is, to some extent, sub¬ 
jective: the “rectangle” in Fig ^ was drawn in order to include 
the already TeV detected sources and sources like Mkn 501 


and Mkn 421 if they were at a redshift ~ 0.1. In Fig. ^ an in¬ 
crease of the redshift corresponds to decrease the fluxes along 
the lines of constant ofRx (the changing K-correction causes a 
negligible deviation, on this scale). 

Fig. ^ shows how the BL Lac objects are placed in the op¬ 
tical - X-ray flux plane. As the seed photons most effective 
for the TeV emission are in the IR-optical band, the optical 
flux could be in principle a better indicator of the density of 
seed photons than the radio flux. However, due to the pos¬ 
sible contamination (either as emission and absorption) from 
the host galaxy (and uncertainties on the intervening reddening 
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Name 

S 

Ex 

Ra 

R 

IJ 

oth. 

z 

Ts GHz 

-P5500A 

ElkeV 

Ref 









Jy 

mjy 

pjy 


0033+595 

X 



X 



0.086“ 

0.066 

2.95 

5.66 

CoOl, SG99, Br97 

0120+340 

X 



X 



0.272 

0.045 

1.23 

4.73 

CoOl, BS94 

0136+391 



X 

X 



— 

0.049 

3.41 

2.39 

BaOO, Br97 

0214+517 



X 

X 



0.049 

0.161 

0.39 

1.60 

Ne94 

0219+428 

X 

X 





0.444 

0.806 

5.12 

1.57 

Pi93, Si91, RoOO, WW90, C97, Di96 

0229+200 

X 






0.139 

0.049 

9.93 

2.88 

AhOO 

0323+022 

X 

X 


X 



0.147 

0.042 

0.81 

4.49 

Fo98, SG99, Be01,Ch99 

0414+009 

X 

X 


X 



0.287 

0.070 

0.82 

9.29 

WoOO, AhOO 

0548-322 

X 

X 





0.069 

0.170 

3.61 

7.47 

CoOl, GAM95 , WW90, Ta95, Ti94, Ro99, Bk99 

0556-384 







0.034 

0.068 

0.56 

2.2 

SG96, GT96, Ci95 

0647+250 

X 



X 



— 

0.073 

5.12 

6.01 

SG99, GAM95, BS94, AhOO 

0806+524 

X 



X 



0.138 

0.172 

6.33 

3.51 

SG99, GAM95, BS94 

0809+024 






X 

— 

0.071 

6.23 

1.85 

BaOO, Br97 

0851+203 

X 

X 


X 

X 


0.306 

1.740 

14.9 

1.70 

To96, Br89, IN88, Ge94, 

E94, Li94, Sa94, Ma88, POl 

1011+496 

X 


X 

X 



0.200 

0.286 

2.04 

1.38 

Fo98 

1028+511 

X 


X 

X 



0.361 

0.044 

1.29 

4.80 

Fo98, BeOl 

1101-232 

X 






0.186 

0.066 

2.52 

9.25 

WoOO, SG99, Ch99 

1114+202 






X 

0.139 

0.074 

3.03 

7.31 

CoOl, Br97 

1133+704 

X 

X 

X 

X 



0.045 

0.274 

9.75 

4.56 

IN88, La96, PirOO, Bk99 

1136.5+6737 



X 

X 



0.135 

0.040 

0.49 

3.40 

F96, Ba94 

1215+303 

X 

X 

X 

X 



0.237 

0.445 

3.08 

3.7 

GAM95, Fo98, Te98 

1218+304 

X 

X 

X 

X 



0.182 

0.056 

1.63 

6.39 

CoOl, GAM95, Pi93, Sa94, Fo98, PirOO 

1417+257 



X 

X 



0.237 

0.040 

0.56 

2.65 

BaOO, La96 

1440+122 

X 



X 



0.162 

0.041 

0.90 

1.47 

BaOO, Br95 

1544+820 

X 






— 

0.043 

1.15 

2.31 

GAM95, Fo98, BeOl 

1553+113 

X 



X 



0.360 

0.636 

6.15 

6.54 

GAM95, Ne94, Fa90, IN88, BeOl 

1722+119 


X 

X 

X 



0.018 

0.088 

2.95 

4.02 

Be92, Pi93, GAM95, SG99, Sa94 

1727+502 

X 

X 


X 



0.055 

0.159 

1.27 

2.73 

Fo98, Pi93, Ke95, Bk99 

1741 + 196 

X 



X 



0.084 

0.223 

1.86 

2.88 

Br97, PirOO 

1959+650 

X 






0.047 

0.252 

1.35 

9.29 

GAM95, SG99, F96, BeOl, Bk99 

2005-489 

X 

X 



X 


0.071 

1.192 

4.84 

5.42 

TaOl, Ki99, Bk99 

2200+420 

X 

X 


X 

X 


0.069 

3.593 

17.27 

1.91 

St94, Br89, IN88, To96, Li94, Bk99, 

Ca97, Sa99, Ma99, POl 

2356-309 






X 

0.165 

0.065 

0.67 

5.78 

CoOl, Be92, Fa94 


“ tentative redshift, Perlman priv. comm, (see NED notes). 


Table 2. TeV BL Lac candidates. S:Slew survey sample; Ex: Exosat sample; R: RGB sample; Ra: RASS sample; IJ: 1 Jy sample; 
oth.: others catalogues, see Donato et al. 2001 and references therein. None of the above sources is in the EMSS sample. The 
references report the basic source of data for the SEDs, in addition to those in the NED database and the respective sample 
papers (Sect. 4.1). All the MeV-GeV upper limits are from Eichtel et al. 1994 (first EGRET catalogue). AhOO: Aharonian et 
al. 2000; Ba94: Bade et al. 1994; BaOO: Bauer et al. 2000; Be92: Bersanelli et al. 1992; BeOl: Beckman et al. 2001; Bk99: 
Buckley 1999; Br89: Brown et al. 1989; Br95: Brinkmann et al. 1995; Br97: Brinkmann et al. 1997; BS94: Brinkmann & Siebert 
1994; C97 : Comastri et al. 1997; Ca97: Catanese et a. 1997a; Ch99: Chadwick et al. 1999; Ci95: Ciliegi et al. 1995; CoOl: 
Costamante et al. 2001; Di96: Dingus et al. 1996; E94: Edelson et al. 1994; E96: Eruscione 1996; Ea90: Ealomo et al. 1990; 
Ea94: Ealomo et al. 1994; Eo98: Eossati et al. 1998; GAM95: Giommi, Ansari & Micol 1995; Ge94: Gear et al. 1994; GS95: 
Ghosh & Soundararajaperumal 1995; GT96: George & Turner 1996; IN88: Inpey & Neugebauer 1988; Ke95: Kerricket al 1995; 
Ki97: Kifune et al. 1997; La96: Lamer et al. 1996; Li94: Litchfield et al. 1994; Ma88: Madejski et al. 1988; Ma99: Madejski 
et al. 1999; Ne94: Neumann et al. 1994; Pi93: Plan et al. 1993; PirOO: E. Piron PhD thesis, 2000; POl: Padovani et al. 2001; 
Ro99: Roberts et al. 1999; RoOl: Robson et al. 2001; Sa94: Sambruna et al. 1994; SG99: Stevens & Gear 1999; St94: Stevens 
et al. 1994; Si91: Sitko & Sitko 1991; TaOl: Tagliafem et al. 2001 and references therein; Te98: Terasranta et al. 1998; To96: 
Tornikoski et al. 1996; WoOO: Wolter et al. 2000 and references therein; WW90: Worrall & Wilkes 1990. 


medium), we consider the radio flux a more reliable indicator 
of the low-energy non-thermal nuclear emission from these 
objects. Note that, in any case, all but one sources within the 
radio-X-ray rectangle are also within the optical-X-ray one. 
Eor a given X-ray flux, then, a relatively high emission in both 


the radio and optical bands should represent a reliable indica¬ 
tion of a large energy density of seed photons. 

Table 2 reports the list of the objects which are inside the 
two “rectangles” in Pig. ^ and Pig. according to the values 
reported in the respective samples (see Section 4.1), and that 
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Name 

L' 

erg 

R 

cm 

B 

G 

r 

e 

n 

71 

O^peak 

72 

rpb 

-^(>40 GeV) 

rpb 

-^(>0.3 TeV) 

rpb 

-^(>1 TeV) 

0033-f595 

1.0e41 

9.0el5 

0.8 

14 

3.0 

3.01 

1000 

5.3e4 

6.0e5 

15.0/2.93 

2.04/0.25 

0.48/0.04 

0120-f340 

4.4e41 

1.0el6 

0.6 

15 

3.0 

3.2 

700 

7.0e4 

2.5e5 

3.17/4.90 

0.28/0.30 

0.06/ — 

0136-f391“ 

3.0e42 

1.0el6 

2.0 

15 

4.8 

3.6 

500 

4.5e3 

5.0e5 

5.22/4.00 

0.56/0.14 

0.12/2.7e-3 

0214-b517 

5.0e40 

8.0el5 

0.4 

13 

5.0 

3.95 

1500 

1.7e5 

6.0e5 

43.4/1.59 

5.93/0.07 

1.43/6.2e-3 

0219-f428 

5.0e42 

2.0el6 

1.4 

15 

3.0 

3.4 

700 

4.7e3 

8.0e4 

7.01/9.62 

0.14/ — 

0.01 / — 

0229-f200 

1.5e41 

1.0el6 

0.3 

15 

3.5 

3.5 

700 

2.8e5 

4.0e5 

7.67/3.81 

0.96/0.31 

0.21 /4.0e-3 

0323-f022 

1.0e41 

1.0el6 

0.9 

11 

3.7 

3.5 

800 

2.7e4 

1.8e5 

6.65/1.18 

0.84/0.01 

0.18/ — 

0414-b009 

8.0e41 

1.0el6 

1.5 

14 

3.5 

3.2 

500 

9.8e3 

1.6e5 

2.91/3.42 

0.23/0.07 

0.04/ — 

0548-322 

7.5e40 

8.0el5 

0.8 

12 

4.0 

3.1 

700 

5.0e4 

5.0e5 

31.9/1.56 

4.14/0.10 

0.91 / 0.015 

0556-383 

7.0e40 

8.0el5 

0.4 

16 

5.0 

3.2 

1500 

2.5e5 

2.5e5 

37.8/5.51 

5.84/0.42 

1.56/ — 

0647-b250“ 

1.0e42 

1.2el6 

1.4 

16 

3.3 

3.6 

500 

1.0e4 

2.0e5 

6.16/8.74 

0.59/0.24 

0.12/ — 

0806-f524 

1.0e42 

1.2el6 

1.5 

15 

4.0 

3.4 

300 

8.5e3 

7.0e4 

14.7 / 10.7 

1.36/ — 

0.27/ — 

0809-f024“ 

8.0e41 

1.0el6 

1.3 

14 

4.8 

3.3 

1000 

1.2e4 

2.0e5 

6.08 / 2.20 

0.58/0.04 

0.12/ — 

08514-202 

5.0e42 

1.3el6 

5.0 

13 

3.5 

3.4 

150 

6.5e2 

7.e3 

23.7/ — 

0.42/ — 

0.03/ — 

1028-4511 

2.7e41 

1.0el6 

1.0 

14 

3.0 

3.2 

800 

2.6e4 

2.0e5 

7.15/3.89 

0.43/ — 

0.06/ — 

1011-4496 

7.0e41 

2.0el6 

0.7 

12 

4.0 

3.4 

300 

1.9e4 

1.0e5 

1.67/3.31 

0.12/0.14 

0.02/ — 

1101-232 

5.3e41 

7.0el5 

0.9 

16 

2.8 

3.1 

300 

5.3e4 

1.0e6 

6.67/10.2 

0.67/0.93 

0.15/0.18 

1114-4202 

4.5e41 

8.0el5 

1.5 

16 

3.0 

4.6 

6000 

6.0e3 

2.5e5 

10.1/8.51 

1.17/0.10 

0.28/ — 

1133-4704 

2.0e41 

2.0el6 

0.8 

10 

5.0 

3.7 

400 

1.6e4 

1.7e5 

62.8/2.15 

8.50/0.03 

1.93/ — 

1136-4673 

4.0e41 

1.0el6 

1.0 

15 

4.2 

3.5 

1000 

2.4e4 

1.5e5 

7.30/5.40 

0.92/0.10 

0.21 / — 

1215-4303 

1.0e41 

1.0el6 

2.5 

11 

4.6 

3.8 

300 

3.4e3 

4.0e4 

4.06 / 0.07 

0.16/ — 

0.02/ — 

1218-4304 

2.0e41 

1.0el6 

1.5 

16 

2.9 

3.9 

600 

1.4e4 

3.0e5 

6.36/5.82 

0.67/0.16 

0.15/ — 

1417-4257 

1.7e42 

1.0el6 

1.5 

15 

4.2 

3.4 

400 

8.4e3 

2.5e5 

3.76/6.88 

0.38/0.21 

0.08/ — 

1440-4122 

2.6e41 

1.2el6 

0.7 

15 

4.2 

3.5 

500 

4.6e4 

5.0e5 

6.11/1.89 

0.78/0.09 

0.20/0.01 

1544-4820“ 

6.0e41 

8.0el5 

1.2 

14 

4.0 

3.3 

400 

1.8e4 

4.0e5 

4.89/5.37 

0.54/0.22 

0.12/ — 

1553-4113 

2.5e42 

3.0el6 

0.7 

15 

2.5 

3.6 

300 

1.3e4 

1.0e5 

8.92/22.3 

0.20/0.42 

0.02/ — 

1722-4119 

4.0e40 

1.5el6 

0.7 

10 

6.5 

3.6 

600 

2.6e4 

5.0e5 

76.6/1.06 

12.8/0.015 

3.52/1.0e-3 

1727-4502 

l.le41 

1.0el6 

0.8 

10 

5.0 

3.5 

600 

3.0e4 

1.9e5 

38.7 / 2.64 

5.19/0.07 

1.23/ — 

1741-4196 

1.0e41 

1.0el6 

0.4 

13 

4.0 

3.8 

400 

1.3e5 

6.0e5 

31.6/4.31 

3.59/0.29 

0.84/0.01 

1959-4650 

8.0e40 

1.0el6 

1.2 

13 

4.0 

3.6 

500 

1.9e4 

1.5e5 

56.7 / 2.08 

7.46/0.03 

1.74/ — 

2005-489 

3.5e41 

1.0el6 

1.2 

16 

3.0 

3.3 

600 

2.1e4 

5.0e5 

67.1 /62.5 

5.14/2.67 

0.90/0.17 

2200-4428 

1.3e42 

5.0el5 

1.2 

14 

3.3 

4.3 

1000 

1.0e3 

2.0e5 

67.7/42.8 

3.32/0.17 

0.43/ — 

2356-309 

2.0e41 

7.0el5 

1.2 

16 

3.0 

3.01 

300 

6.0e4 

7.0e5 

7.64/3.30 

0.84/0.19 

0.12/0.03 


Table 3. Input parameters for the SSC model, resulting values of 7 peak, and predicted fluxes at high energies, above 40 GeV, 0.3 
TeV and 1 TeV, according to the parameterization of the SED adapted from Fossati et al. (1998) (first number) and according 
to the SSC model discussed in Section 5.1 (second number). L' is the intrinsic power (i.e. measured in the comoving frame), R 
the cross sectional radius of the emitting region, T is the bulk Lorentz factor, 9 is the viewing angle, n the slope of the particle 
distribution above the cooling energy 7 c (see text), 71 and 72 are the extreme Lorentz factors of the injected particle distribution 
and 7 peak is the particle Lorentz factor of the electrons emitting most of the radiation (i.e. at the synchrotron and self-Compton 
peaks). z = 0.2 assumed. ^ Fluxes in units of 10“^^ photons cm“^ s“t. 


therefore we consider the best candidates for a possible TeV 
detection. 


5. Prediction of the TeV flux 

To better quantify the predicted high energy flux of our best 
candidates we have collected from the literature other data for 
all the objects listed in Table 2, in order to construct their SED. 
We have then used two different methods to estimate their high 
energy emission. First, for each source, we have applied a SSC 
model, as explained below, aimed to fit the synchrotron compo¬ 
nent of their SED and to predict the inverse Compton spectrum. 
Then we have also calculated the predicted spectrum according 
to a slightly modified version of the parameterization given by 
F98, thought to describe the average SED of blazars. 


Note that the predicted SEDs shown in Fig. 6 and the TeV 
fluxes listed in Table 3 do not take into account the absorption 
of TeV photons by the IRB, since this is indeed one of the im¬ 
portant unknowns we would like to measure. 

5.1. Homogeneous Synchrotron Self-Compton 
model 

We have applied a homogeneous, one-zone synchrotron self- 
Compton model to our best TeV candidate BL Lacs. This 
model is “one-zone version” of the model in Spada et al. 
(2001) and is described in detail in Ghisellini, Celotti & 
Costamante (2001). The main characterizing feature of this 
model is the assumption of an particle injection mechanism 
lasting for a finite time. In this case the emitting particle dis¬ 
tribution never reaches a complete steady state. A physical see- 
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nario where such a behavior naturally occurs is that of inter¬ 
nal shocks (see e.g. Piran 1999; Ghisellini 1999; Spada et al. 
2001 ), i.e. collisions of different parts of the jet plasma moving 
at slightly different speeds, naturally leading to shocks lasting 
for a finite time. This scenario is suggested by the rapid vari¬ 
ability always present in BL Lac objects in general and HBL in 
particular, especially at high energies. We briefly outline here 
the other main assumptions of the model. 

We assume that the jet is conical, with half opening angle 
~ 1 /r, and approximate the emitting region as a cylinder, of 
radius R (= 'ijjz, where z is the distance along the jet axis) and 
width Ai?' = R/T (in the comoving frame, here L is the bulk 
Lorentz factor). This corresponds to assume Ai?' constant up 
to the first collision, in the internal shocks scenario (Spada et 
al. 2001). We derive the particle distribution assuming that the 
acceleration mechanism is in the form of a continuous injection 
of relativistic particles distributed in energy as a power law of 
index n—1 between 71 and 72 . The injection is assumed to last 
for a finite time, set equal to tinj = AR'/c. We define 7 c as the 
energy of those electrons that can cool in the injection time fjnj, 
i.e. fcooi( 7 c) = Lnj-At energies greater than 7 c, particles radia- 
tively cool, and the distribution reaches a steady state in a time 
smaller than finj. As a consequence, the emitting particle distri¬ 
bution is assumed to be a power law of index n above 7 c. Below 
this value, there can be different cases according if 7 c is greater 
or smaller than 71 . If 7 c > 71 , we have N{j) oc 7 “("“i) be¬ 
tween 7 i and 7 c. Alternatively, if jc < 7i, then JV(j) ex 7 “^ 
between 7 c and 71 . We further assume that, below the mini¬ 
mum between 71 and 7 c, ^( 7 ) ex 7 “^. According to these 
assumptions, the random Lorentz factor 7 peak of the electrons 
emitting most of the radiation (i.e. emitting at the peaks of the 
SEDs) is determined by the importance of radiative losses and 
can have values within the range 71 - 72 . Its value is listed in 
Table 3. The source is assumed to emit an intrinsic luminosity 
L' and is assumed to be observed with the viewing angle 0. All 
these input parameters are listed in Table 3, together with the 
predicted high energy photon fluxes above three representative 
frequencies. The resulting fits are shown in Fig. 6 as solid lines. 

In conventional SSC models (Tavecchio, Maraschi & 
Ghisellini 1998), the knowledge of the location of the two 
peaks and their fluxes, together with the variability timescale, 
suffices to derive unambiguously all relevant physical parame¬ 
ters. Here, instead, we lack two important observables, namely 
the frequency and flux of the Compton peak. We therefore need 
to supply two additional relations with respect to the “standard” 
SSC model. 

One of these relations comes from our assumption of finite 
injection time, resulting in a relation between r'peak and the 
magnetic field. In fact, if synchrotron losses are relevant (as 
in the sources already detected at TeV energies, for which the 
high energy component never exceeds the synchrotron one), we 
have a constraint on the value of the magnetic field, controlling 
the radiative losses timescales and hence the value of 7 peak- 

For the second unknown to be provided, we decided to limit 
the value of the beaming Doppler factor S within the relatively 
narrow range 9 < i5 < 20. 

Note, however, that also within the previous assumptions, 
the 7 -ray flux predictions still have a large uncertainty, mainly 


due to the dependence of the flux from the emission volume 
(i.e. from R and F). In our model, however, these are not com¬ 
pletely free parameters, since on one hand they must be com¬ 
patible with the variability timescales, and on the other they 
control the injection timescale finj. The dependence of the in¬ 
verse Compton flux from R and F is therefore complex. For 
instance, for a given L' and B, an increase of R makes the 
Compton flux to decrease (ex R~^), but it makes also tinj to in¬ 
crease, leading to a smaller value of 7 c (= 7 peak), and then to 
a smaller synchrotron peak frequency. To compensate for that, 
it would be necessary to decrease B (to increase the cooling 
time), but this leads to an increase of the Compton flux. 

Our choice has been to use values of R around 10^® cm, 
which are the typical values obtained for sources with good 
monitoring and with SEDs well sampled also at high energies, 
like Mkn421 and lES 2155-304. 

Another point of uncertainty is the exact determination of 
:^peak and the corresponding synchrotron peak flux. For many 
sources this is provided by the BeppoSAX observations, which 
can be fitted by a broken power law, immediately yielding 
v-peak- For Other sources of poorly known SED the determina¬ 
tion of r-peak is more uncertain, and for these sources the pre¬ 
dicted SEDs of course suffer from this uncertainty (see e.g. the 
SED of 0214H-517, 1440H-122 and 1544H-820, for which there 
is no information of the slope of the X-ray spectra and only a 
few observations in the IR-optical bands). 

The applied model is aimed to reproduce the spectrum orig¬ 
inating in a limited part of the jet, thought to be responsible of 
most of the emission. This region is necessarily compact, since 
it must account for the fast variability shown by all blazars, es¬ 
pecially at high frequencies. Therefore the radio emission from 
this compact regions is strongly self-absorbed, and the model 
cannot account for the observed radio flux. 

5.2. Fossati et al. (1998) description of the SED 

Fossati et al. (1998) proposed a simple phenomenological de¬ 
scription of the average SED of blazars based on their bolomet- 
ric observed luminosity, thought to be well traced by the radio 
luminosity. In their parameterization, the radio luminosity is as¬ 
sumed to be the key parameter, determining the peak frequency 
of the synchrotron spectrum and the relative importance of the 
inverse Compton power. More recently, Donato et al. (2001) 
have revisited this parameterization, assuming a slightly dif¬ 
ferent relation between the radio power, the synchrotron peak 
frequency and the inverse Compton luminosity, but only for ob¬ 
jects below a radio luminosity of 10"^^ erg s“^. With this new 
parameterization, objects of low power are assumed to have 
equal luminosities in the synchrotron and in the self-Compton 
components of their spectra, and the ratio of the peak frequen¬ 
cies (t'e/^'s) was allowed to increase for increasing Vg (as ex¬ 
pected in a SSC scenario). 

To improve this parameterization also for objects with even 
lower radio luminosity (i.e. higher r'peak), we have further mod¬ 
ified the Donato et al. prescription for sources with radio lu¬ 
minosity below erg s“^, to mimic some effects of the 

Klein-Nishina regime. The first modification involves the ra- 
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Log V [Hz] 

Fig. 6 . a): SEDs of our best candidates for TeV emission. The solid lines refer to the SSC model as explained in Section 5.1. 
Dashed lines correspond to the phenomenological prescription of Fossati et al. (1998), as slightly modified by Donato et al. 
(2001) and in this paper. Sources of data listed in Table 2. 


tio Vcjvs, since in the Klein-Nishina regime, Vg oc 7pea,k 
while Vc oc 7peak- Then the ratio v^/vs oc 7 p/g^,^ oc . 

Below Lr = 10^^'^ erg s“^, we then decrease Vc/vs assum¬ 
ing that \og{vc/vs) = 9.4 — 0.8 x [41.2 — logLR]. The other 
modihcation concerns the width of the parabola representing 
the Compton peak, which is reduced with respect to the syn¬ 
chrotron one (cTcomp = '^sync/2; See Fossati et al. (1997), ap¬ 
pendix A, and F98 for a description of the parameters). In this 
way the resulting SED shapes provide a better representation 
of the individual SEDs for the highest t'peak sources. 

We have applied this parameterization to all the objects in 
our sample, and the resulting hts are shown as dashed lines in 


Fig. 6 . From these we obtained the photon fluxes listed in Table 
3, integrating above the corresponding threshold energy. 

From the comparison with the homogeneous SSC models, 
we can see that the F98 parameterization tends to overestimate 
the Compton 7 -ray emission, although it agrees with the ex¬ 
isting upper limits in the GeV and TeV bands in all cases but 
0851H-202, 1133H-704 and 1553H-113. This is not unexpected, 
since, by construction, the Compton power is never less than 
the synchrotron one and the assumed parabolic shape does 
not account for sharp high energy cut-offs. This parameteri¬ 
zation however, which is built to describe the average SED of 
sources of equal synchrotron and self-Compton power, is in 
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Fig. 6. b); SEDs of our best candidates for TeV emission. The solid lines refer to the SSC model as explained in Section 5.1. 
Dashed lines correspond to the phenomenological prescription of Fossati et al. (1998), as slightly modified by Donato et al. 
(2001) and in this paper. Sources of data listed in Table 2. 


good agreement with the flux measurements of the already de¬ 
tected TeV sources (see Fig. |^. 

6. Discussion and conciusions 

With new Cherenkov telescopes foreseen to operate in the next 
few years the TeV extragalactic astronomy is entering its adult¬ 
hood. A tenfold increase in sensitivity (expected for the forth¬ 
coming installations) would mean the possible detection of ~ 
100 BL Lacs, if the counts of BL Lac object at TeV energies 
are roughly Euclidean (in the bright flux end) and neglecting 
absorption by IRB. Because of IRB absorption, the counts will 
be flatter than Euclidean, but the lower energy threshold of 
some new instruments may compensate for the extragalactic 


absorption, as well as favoring the detection of slightly less 
blue objects. Therefore many more sources are expected to be 
detectable by the new telescopes, and this motivated us to study 
which kind of BL Lac objects is more likely to be detected at 
high energies. 

Our findings can be summarized very simply once we 
realize that, to produce a strong TeV emission, the inverse 
Compton process needs a sufficient number of both very high 
energy electrons and soft seed photons. Therefore we require 
both a strong X-ray flux and a sufficiently strong radio- 
through-optical flux. Since the optical flux can be contami¬ 
nated, especially in low redshift sources, by the underlying host 
galaxy, our sources are primarily selected as bright both in the 
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Fig. 6. c): SEDs of our best candidates for TeV emission. The solid lines refer to the SSC model as explained in Section 5.1. 
Dashed lines correspond to the phenomenological prescription of Fossati et al. (1998), as slightly modified by Donato et al. 
(2001) and in this paper. Sources of data listed in Table 2. 

X-ray and radio bands. All but one of these sources are also 
bright in the optical. 

With respect to the previous work by Stecker et al. (1996), 
our criterium introduces the further requirement that the source 
must be a relatively strong radio emitter. The other difference 
is that we considered not only the Einstein Slew survey sample 
of BL Lacs, but several other BL Lac samples. 

Besides selecting the best candidates through their location 
in the radio - X-ray flux plane (a criterium largely model in¬ 
dependent), we have also tried to quantify the level of the ex¬ 
pected high energy emission for each selected source, by ap¬ 
plying a one-zone synchrotron self-Compton model and also 
the phenomenological description of the SED of Fossati et al. 


(1998), slightly modihed to better account of the average SED 
of low power BL Lacs. The latter model, by construction, as¬ 
sumes equal power between the synchrotron and the inverse 
Compton components of the SEDs, and almost always pre¬ 
dicts larger high energy fluxes than the SSC model. In the 
SSC model, in fact, the Compton dominance is not fixed a 
priori, but found by fitting the synchrotron part of the spec¬ 
trum, which Axes the value of the magnetic held. We stress that 
the F98 prescription was designed including also the SEDs of 
those BL Lacs already detected in the TeV band (and indeed 
is in good agreement with their TeV flux levels), and therefore 
seems more appropriate to predict the TeV flux of sources in 
high state. The adopted SSC model, instead, is designed to fit 
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Fig. 6. d); SEDs of our best candidates for TeV emission. The solid lines refer to the SSC model as explained in Section 5.1. 
Dashed lines correspond to the phenomenological prescription of Fossati et al. (1998), as slightly modified by Donato et al. 
(2001) and in this paper. Sources of data listed in Table 2. 


the known synchrotron part of the SED, which is often repre¬ 
sentative of a more “normal” or quiescent state. This explains 
the sometimes large discrepancy between the predicted fluxes 
of the two models. Since BL Lac objects are among the most 
variable sources, especially at high energies, the two foreseen 
flux levels could be thought of as an approximate range of vari¬ 
ability, and the average flux could be considered as a measure 
of the probability to And the source in a particular TeV state. 

We would like to stress, anyway, that the uncertainties 
on the key parameters we used for the model and the non¬ 
simultaneity of the fitted data lead to large uncertainties in the 
predicted TeV flux, sometimes of the same order of the differ¬ 
ences between the two adopted models. Our predicted fluxes, 
therefore, also in the case of the SSC model, must be consid¬ 
ered as “best guesses” on the high energy emission from these 
objects. Within these limits, the SSC model provides more in¬ 
formation than the phenomenological parameterization, since 
it gives also the expected shape of the high energy spectrum. 

Since the level of the synchrotron X-ray flux measures, in 
our scenario (as well as in the Stecker et al. 1996 one), the 
number of TeV energy electrons, the X-ray monitoring of our 
candidates is particularly useful to catch sources in high TeV 
states, as already partially done through the All Sky Monitor 
(ASM) onboard the RossiXTE satellite. 


Besides the fluxes above 300 GeV, the most common 
threshold of present Cherenkov telescopes, we also give our 
estimates above 40 GeV, which is approximately the energy 
threshold of CELESTE and of forthcoming observatories like 
HESS and MAGIC (with VERITAS at ~ 50 GeV, Weekes 
1999). Emission at these energies is much less absorbed by 
the cosmic infrared background, giving the opportunity to see 
more distant sources and study their intrinsic spectrum in an 
unabsorbed band. In addition this energy range will link ground 
based Cherenkov observations and the data coming from satel¬ 
lites, such as AGILE and GLAST, observing from a few tens 
of MeV to a few tens of GeV. 

As a final note, we warn that our flux estimates do not in¬ 
clude the possible absorption due to the infrared background, 
since we preferred to be independent of this factor. In fact, we 
focussed on the conditions for the TeV emission, providing a 
list of possible sources, in order to allow an independent test of 
the IR absorption effects. This differs with respect to the flux 
estimates in Stecker et al. (1996), which instead account for the 
IR background absorption. However, since most of the photon 
flux is expected to be emitted at energies below 1 TeV (as can 
be seen in Table 3, comparing the values above 0.3 and 1 TeV), 
the reported fluxes should not be much affected for sources up 
to z ~ 0.1, according to the present estimates on the IR back¬ 
ground (see Stecker 2001 and references therein). 
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